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Abstract 

The gyrokinetics formulation of plasmas in strong magnetic fields 
aims at the elimination of the angle associated with the Larmor rota- 
tion of charged particles around the magnetic field lines. In a pertur- 
bative treatment or as a time-averaging procedure, gyrokinetics is in 
general an approximation to the true dynamics. Here we discuss the 
conditions under which gyrokinetics is either an approximation or an 
exact operation in the framework of reduction of dynamical systems 
with symmetry. 



1 Introduction 

In contrast with the many types of approximation schemes used to deal with 
complex problems, a reduction is an exact process. Reduction is a process 
by which, given a multidimensional dynamical system, one focus on a small 
subset of variables, obtaining exact equations for the variables in this subset. 
Along the way, of course, one looses information about some of the dynamical 
details of the system. In particular, several distinct evolutions of the whole 
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system may lead to the same dynamical trajectory when projected on the 
space of the reduced variables but, in any case, one is confident that the 
dynamics of the reduced variables is exact. 

Reduction in mechanics had its origins in the classical works of Euler, La- 
grange, Hamilton, Jacobi, Routh and Poincare. Routh's elimination of cyclic 
variables and Jacobi's elimination of the node are among the first examples 
(for reviews of the historical aspects and the several types of reduction we 
refer to [I] [2] In general, reduction applies when the system under 

study possesses some kind of symmetry. Then, the variables associated to 
the invariance directions are eliminated and one ends up with a dynamical 
description on a quotient space. 

In the context of symplectic manifolds, a construction has now become 
standard [I] [3] which, for reference and later use, we now summarize: 

Let (P, uj) be a symplectic manifold where a Lie group G acts by symplec- 
tomorphisms (<j>gUJ = u) . An equivariant moment map for the group action 
is a map J : P —>■ g* (g* being the dual of the Lie algebra g of G) such that 

if J denotes the dual map from g to the space of smooth functions on P, we 
have 

d \J (0 1 = Hp" (1) 

and J ((fig (x)) = Ad*-! (J (x)) , Vx G P, g G G. If the momentum map is not 
equivariant it can be converted into an equivariant one by central extension 
of the group |5J. Eq.([T]) means that each infinitesimal generator £p of g has 

J (£) as an Hamiltonian function. 

Let /i be an element of g* and G^ its coadjoint isotropy subgroup. For 
any G— invariant Hamiltonian, J -1 (/i) C P is an invariant set for the dy- 
namics and the reduced space J -1 (fi) /G^ = P M is a symplectic manifold 
with symplectic form f2 M determined by 7r*^ = tt^ is the projection 

J -1 (/i) — > P^ and the inclusion J -1 (/z) — > P. 

In the symplectic reduction the emphasis is on the projection of the Pois- 
son structure and Hamiltonian dynamics to a quotient space by the action of 
a symmetry group. In the Lagrangian approach, reduction follows a different 
approach. Reduced variables are identified and then one proceeds to carry 
the variational structure to a quotient space. In the process not only the com- 
pleteness of the reduced variables set should be checked, but also whether 
the variational structure can be carried over to the quotient space. Here, in 
principle, a symmetry group need not to be known to begin with. One may 
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start by conjecturing some tentative set of reduced variables (functions of 
the original variables). Then, using the Poisson brackets of the original vari- 
ables, check whether this trial set is algebraically closed. If not, the Poisson 
brackets themselves will suggest new variables to close the set. Then, to be 
sure that the variational structure is carried over to the quotient space, it 
is necessary to check whether the (exact) dynamics of the reduced variables 
may be obtained from a Lagrangian written on these variables alone. If that 
is possible, an exact reduction has been performed. 

Even when an exact reduction on some subset of variables is not pos- 
sible, the reduction point of view provides a step by step approach to an 
approximation scheme. For example, it is possible that the trial set of re- 
duced variables is algebraically closed with a consistent bracket structure, 
but that their dynamics cannot be entirely defined by an action principle 
containing only these variables. To proceed it must be necessary at this step 
to propose some kind of approximation, but at least the method makes it 
clear how much of the reduction process is exact and how much it implies an 
approximation. 

The same methodology of approximation through reduction may be used 
in the group theoretical setting. Suppose that some dynamical process is 
known that is not exactly symmetric, but we are interested only on the 
dynamical features consistent with the symmetry. Then an approximation- 
through-reduction may be obtained by projecting the dynamics on the space 
of functions that possess the required symmetry. 

In the context of plasma physics in strong magnetic fields, two very dif- 
ferent time scales rule the physical phenomena. One is the fast time scale 
of gyromotion of the ions around the magnetic field (of order 27r ^ c ) and the 
other the longer time scale associated to the electric and magnetic gradient 
and curvature drifts. An important simplification arises when the two time 
scales can be treated separately. This might not be possible, for example in 
the study of fluctuation phenomena. However, for some stability and trans- 
port problems it is indeed an useful approximation. When the fast time 
scale is separated, or integrated over, one obtains the so called gyrokinetic 
and gyrocenter equations. The usual way this is done is either by simple 
averaging the single particle dynamics over the gyroperiod [6] [7] or, more 
accurately, by performing a Lie-transform perturbative expansion to obtain 
an action two-form where the gyroangle dependence is to be asymptotically 
eliminated. In both cases the point of view is to address the gyrokinetics for- 
mulation as an approximation to the plasma dynamics. The Lie-transform 
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perturbative approach has been extensively developed [S] [S] [ID] [H] leading, 
for example, to a third-order perturbative analysis of a plasma moving with 
a nonuniform fluid velocity [12] . 

Here, our purpose is to find out how much of the gyrokinetics program 
may be framed as an exact reduction and how much it is an approximation. 
In Sect. 2 we deal with Lagrangian reduction as a time-averaging process, fol- 
lowing a procedure quite analogous to the one that leads from the Lagrangian 
dynamics of many particles to the Euler equation. The only difference is that 
now we define time-averaged Eulerian coordinates. 

The Lie-transform approach to gyrokinetics ([13J and references therein) 
has been extensively studied and provides successive approximations of prac- 
tical interest to the reduced Vlasov and Maxwell equations. Its purpose is to 
obtain the asymptotic elimination of the gyroangle dependence. However to 
show that this procedure provides an exact reduction of the dynamics would 
require the proof of convergence of the perturbative series, what has not yet 
been done. In Sect. 3, following a non-perturbative approach, we show how 
gyrokinetics may be considered an exact reduction of the dynamics in the 
sense of the symplectic reduction scheme discussed above [I]. 

2 Lagrangian time-averaging reduction 

Here we discuss, gyrokinetics reduction as a time-averaging process, in partic- 
ular to exhibit the meaning and shortcomings of such an approach. We start 
from the non-relativistic particles plus electromagnetic fields Lagrangian [1] 



L = Jdz f(z ) j f (z , t)-e J dx(^> (x, t) - T ^ * A (x, t) j 5 (x - r (z , *)) 

+ _L J dx(\E(x,t)\ 2 -\B(x,t)\ 2 ) (2) 

zo = ( r o,Po) stands for a particle coordinate in phase space, r(zo,t) is the 
position of the particle that at time zero was at zq, and / (zq) is the particle 
density at time zero, that is 

JV 

f(zo) = J2 6 ( r °- r o ) ) S (p°-Po ) ) (3) 

i=l 
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for iV particles. We consider only one ion species with mass m and charge 
e. Generalization of all statements to an ensemble of different ion species is 
straightforward. 

The Lagrangian (jSJ) contains information about the identity and dynamics 
of each particle at all times. We now define our reduced variables 

f h (R,P,t) = J dz f(z )8(R-r(z ,t))d(p-^^) (4) 



where 



r(z ,t) = — I dr)ht(r))r(z ,t-r)) (5) 



ir(z ,t) = — / dr)h t (r))Tr(z ,t-r)) 



A 



The kernel ht (if) is a function equal to one in the interval [— y, y] and zero 
otherwise. The reason for the division by / (zq) in the second delta function 
of Eq.(jH) comes from the fact that the canonical momentum is 

k (*b, t) = — = / (zo) (m r (z , t) + -A (r (z , t))) (6) 

8r(z ,t) V c J 

In the new variables two types of information are lost. First, the parti- 
cle identities are lost, being replaced by densities at points in phase space. 
Second, the dynamics itself is averaged around each time t with a window of 
size At, that need not be the same for all t. 

The next step is to check the consistency of this reduction. First one 
checks whether the new set {fh (R, P, t)} is algebraically closed and whether 
its Poisson bracket may be expressed on the variables (R,P). We denote 
Z = (R, P) and compute the Poisson bracket of the / s in the Lagrangian 
variables (r, 7r). 

{fh(Z,t)f h (z',t)} {r ^ = 

= ^J dz f (zo) [5' (f (z , t )-R)s{j^-P)s(r- R') 5' - P') 

-5 (F (zo, t) - R) 5' - P) 5' (f - K) 5 (jfy - P'X 
= £ [f h (R\ P, t) 5' (R' - R) 5' (P - P') - f h (R, P', t) 5' (R - K) 5' (P' - P)] 
= ijdZ"f h (Z",t){5(Z-Z"),5(Z'-Z")} {Rp) 

(7) 
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One sees that not only is the set {fh (R,P,t)} algebraically closed, but also 
its Poisson bracket may be written in the variables (R,P). The conclusion 
is that this set of (Eulerian) time-averaged densities is a kinematically con- 
sistent reduction. The quantities R and P are in some sense coordinates in 
a time-averaged phase-space. However, one should remember that, from the 
reduction point of view it is the / s that are the actual dynamical variables 
which, with its bracket and reduced Hamiltonian, control the dynamics. The 
set (R, P) is just a set of labels for the dynamical variables. 

The next step in the reduction process is to check the dynamical law and, 
in particular, whether the dynamics may be completely expressed in the new 
variables. The Hamiltonian that follows from ([2]) by the Legendre transform 
is 



H 



/w(*o){^ 



- -A(r {z ,t)) 



f (zb) c 
+i- J dx(\E(x,t)\ 2 + \B(x,t)\ 2 ) 




Then we compute the Poisson bracket in the original Lagrangian coordinates. 
Averaging the computation over the interval r£ [t — tt? t + 4f] > considering 
at each r a complete set of canonical variables (r (zq, t) , it (zq, t)). 

f h (Z, t) = {f h (Z, t) , H} (r ^ = ijdr {f h (Z, t) , H (t)} { , tU(z0!t)) 
= ^JdrJ dz f (z ) [6' (f (z , t)-R)8(*$-P).± - \A (z , r)) 

+ 5 (r (z , t) - R) 5' (^gf - P) ■ (*g£ - *A (z , r)) • VA - eV0 

(9) 

We see that, by taking averages over r and using the multiplication by the 
delta function in the first term, one obtains the evolution equation in the 
renormalized time 4- 



dh(Z,t) ] fP_ _ ^\ df h (Z,t) ] /_e_ _ e \ ^ _ \ df h (Z,t) 

q j \m mc / dR \mc \ f (z ) c J J dP 

(10) 

In this equation the coefficient of the ^ derivative takes the usual form for 
a Vlasov equation in time- averaged variables but, unless the average of the 
products coincide with the product of the averages, that is not the case for 
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the coefficient of the ^ term. For the reduction to be exact one should be 
able to express the dynamics of the reduced densities fa (Z, t) purely in terms 
of the variables (R, P) and the time-averaged field potentials (0, A). It is at 
this step that the reduction may become an approximation. However, if it 
is possible, by an adequate choice of the time-dependent averaging interval 
At, to insure that the average of the products in ([TO]) may be written in the 
averaged variables, then one has an exact reduction. Otherwise the reduced 
dynamics implies an approximation. 

As an example of a particular situation where time-averaging provides 
an exact reduction consider constant electromagnetic fields with the uniform 
magnetic field B oriented along the z— axis. The particle equations of motion 



arc 



, dv x 



<T = eE x + lv y B (11) 



= eE y — zVr.B 



c 



Then 



v x = § (1 - cos (cut)) Ey + -§ {cE x + v y (0) B) sin (W) 
+v x (0) cos {ojt) 

Vy = --|(1 - cos (ut))E x + j; (cE y -v x (0) B) sin (cot) (12) 

+v y (0) cos (cut) 
v z = ME z + v z (0) 

cu = Choosing A = ^j, the corresponding averaged dynamics is 

— 

v., = -c E " 



y ^_b 
v z = ^tE z + v z (0) 





(fir 
dt 



(13) 



dt 





IE, 



dt m 

and the reduced Vlasov equation 

df df df df e- df 

dt x dx y dy z dz m z dv z 

As expected in an exact reduction, the angle, argument of the trigonomet- 
ric functions, being an ignorable coordinate for the averaged dynamics, the 
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number of relevant coordinates is reduced from six, (x, y, z, v x , v y , v z ), to four, 
(x, y, z,v g ). 

Here we have analyzed reduction as a time-averaging process to empha- 
size its different aspects, namely the structural (or kinematical) level that 
concerns the Poisson brackets and the dynamical level that refers to the 
equations of motion. Although we have considered an averaging interval X t 
that depends on time, it must be pointed out that, for general magnetic field 
configurations, the gyrokinetics approximation is more general in the sense 
that the averaging depends not only in time but also on the position of the 
particles. In the present Lagrangian reduction scheme it would correspond 
to having a kernel h = h (zq, t). 

For general electromagnetic fields, the difficulties with the time-averaging 
point of view lead naturally to the more appropriate notion of gyroangle in- 
dependence, which is also the point of view of the Lie-transform perturbative 
approach [13]. Here we also follow the gyroangle independence point of view 
but, instead of a perturbative approach, we prove the existence of an exact 
invariant from which an exact reduction is shown to follow. 



3 Gyrokinetics as an exact reduction 

In the extended phase-space (^x, P,t,h^j, the Hamiltonian of a particle mov- 
ing in an electromagnetic field is 



H[x,P,t,h) = ±-(p-- c A) + e<$>-h (15) 



the non- vanishing elements of the Poisson tensor being 

{x\p>} = S ij 
{t,h} = -1 

Changing coordinates to (x,v,t, k^j with 

v = ±(p-ZA 

m \ c 

k = h - e$ 



(16) 



(17) 



leads to ^ 

H (x,v,t,kj = -m v —k (18) 
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and 

I ' J m 
rr tk = S+ h\ = — 1 

(19) 





m 


o tk = {t, k} = 


-1 


a vivj = {v\v j } = 




a v*k = k y = 


m 



B*J = e ijk B k 

It is convenient to decompose the velocity into magnetic field adapted 
components, 

nil = Vi hi 

(20) 

v± = v - b v ■ b 



b= t-^t , for which the equations of motion are obtained from (11811 and (1191) 

B 

by f = {F,H} 



(21) 

f t v ± = ^i+^(«ixB]-ti'P.V](6ll 



For bounded electromagnetic fields with bounded derivatives, the right-hand 
side of the system of equations ( 12T|) is locally lipschitzian. This insures 
existence and uniqueness of the solution for a time interval. 

To prove the existence of an exact reduction of a dynamical system one 
has to identify a symmetry group or, equivalently, the existence of one or more 
invariants. Our method depends on the construction of a formal invariant. 
Existence of such invariants, to all orders in perturbation theory, for charged 
particle motion in a strong magnetic field was first pointed out by Kruskal[T7]. 
Here we attempt an explicit construction of an exact invariant. The technique 
hinges on transforming the equation for the transverse velocity to the form 

^~v± — u±*J = k ^v± — u±j x B — 6 T fx, v,tj— a (ic, v,t\ ^v± — Uj} 

(22) 



dt 

When this is achieved it is easy to see that 



M 



2 

v± - u± 



(23) 
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is a constant of motion provided 

F = exp(-2/ a[x,v,T)dT] (24) 



The results are summarized in the following proposition: 

In the domain of existence of bounded solutions of the system WIS) and 
for bounded and sufficiently smooth electromagnetic fields, there are func- 
tions T (^x , v ,tj , a (x , v ,tj and a transversal vector function (^x,v,tj 
such that the equation [2^1 holds. For fields that have no explicit time de- 
pendence and sufficiently large magnetic field we also have u±_=u± ( x, v 

a — a(x,v) andT = T\x,v 



Proof : Without explicit time- independence of u±, a and T, the result is 
fairly trivial. It suffices to rewrite ( 1221) as 

u± —k u±_ x B +a (x,^v,tj u ± = + a (^x , v ,t^j — k B v ± + bT (^x , v ,t 

(25) 

and then, the assumed boundedness of the fields and the solutions of (!2~TT) 
implies the existence of a solution u± (^x, v , t^j for each choice of T (^x, v , t^j 

and a (x,v,tj. 

Less trivial is to show the existence of time-independent solutions. Be- 
cause we are going to construct the quantities u±, a and T and the invariant 
as a formal operator series, it is important to control the magnitude of the 
operator action in the space of velocities, in particular to guarantee that it 
does not grow with \B\, the magnetic field intensity. In the second equation 

in (j2T|) we bring the term ^ yu j_ x B^j to the left-hand side 



i + £Bx) v ± =^E±-v t (vV) (bb) (20) 

and rewrite the operator + ^ B x\ as a derivation 

d e\B\ I (i) d (o) d , 
dt mc y ± dvf dv<P 
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t^j 1 and v± being the components of v± in an arbitrary transverse coordinate 
system. Then the equations (1211) become 

Dv\\ = %E\\+ v ± ■ (v -v) b 

V (28) 



We see that the right-hand side no longer involves explicit dependence on 
the , magnetic field, that is, the action of the D operator does not introduce 
frequencies comparable to the Larmor frequency. It is the choice of the 
D operator for the operator expansions that, for large B, implements the 
separation of time scales. 

We now apply the identity, 

aj_= — ^-j (a± xlPjxB (29) 

that holds for transversal fields, to the first and the transversal part of the 
second term in the right-hand side of the second equation in (1281 . to obtain 

D (v ± -c^ + ^(v.v) 6xi)=- (-c^ + ^ (v-v) b x b) x B 

-D {ci^ -^(-v)5xb}- 6 ^ {v -V) b 

(30) 

Separating the transversal and longitudinal components of the term D {■ ■ ■} 
in the right-hand side of Eq. (l3"0l . using the identity (I2"9~j) for the transversal 
component and iterating the process one finally obtains 

d_ 

dt 

with u±, T and a given by 

u ± = fl-#ixD} _1 (c^-5(,.V)6x B ] (32) 



(v± - u±) = + 7) (uj. - w±) x 5 - b r - a (vj. - w ± ) (31) 



and 



+ b -x + v* {v ■ V) b 



(33) 
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x ■ 



and 



a 

X 
( 



u 



± 



u 



7 = - 



dt 



v± - u± 



u 



\ 



B 

\BY 



■ X 



(34) 



(35) 



(36) 



/ 



Notice that 



dt 



dt 



+ [v -V). When D is applied to the velocities it is 



understood that it is a replacement by the right-hand side of Eqs. ff28|) . When 
D is applied to the fields, explicit time dependence in «i and V occurs only if 
the fields themselves have an explicit time dependence. Therefore, whenever 
a proper meaning is given to the series in fl32|) and the fields are not explicitly 
dependent on time, one proves the existence of time-independent solutions 

w_l=m_l yx, vj, T = T vj and a (lie, vj. 

We now discuss the convergence of the formal series in (|32l) . The operator 



B xD having the real line as spectrum, 1- 



e|S|- 



B xD cannot have an inverse 



defined in the whole L 2 . Instead we consider the action of the operator at 
sucessively higher orders on a space of velocities and fields bounded by some 
quantity M. Because of the choice of the operator D, the terms do not grow 
with B, however, because of the nonlinear nature of the ( 1281) action, there is 
a proliferation of terms and, at most, we obtain a bound 



u 



< 



M 

W\ 



which does not insure convergence. However, multiplying and dividing by n\ 
and exchanging the order of sum and integral one obtains 



< 



E 



MzV 
W\) 



where we have used the integral representation 



nl 



dzz n e~ 
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d denoting the integration along the line (1 +s)r] (rj real G [0, oo)) in the 
complex plane. Then 



u 



< 



meaning that u±_ is bounded by a series that is Borel summable along a 

direction not containing the real half-line. Then the formal series in ( I32jl is 
also expected to be Borel summable. In this case, by a theorem of Borel 
[T8] |19j, one solution u exists for which fl32|) is an asymptotic series. The 

" _L 

nature of (|32|) as an asymptotic series has the implication that for practical 
calculations the number of terms to be kept depends on both the magnitudes 
of the velocities and the magnetic field. 

Now from the existence of a vector function u±, satisfying fl22|) . it follows 
that for M, constructed in ( l23l) . = 0. Therefore, an exact reduction 
is now possible following the Marsden-Weinstein theory [I]. M itself is the 

(dual) moment map (the J (£) in EqJTJ that by 

f = W.M} (37) 



generates the action of the symmetry group on phase-space functions Q yx , v 

For each value fi of the invariant M one obtains a symplectic reduced 

space of dimension four. Existence of a set (q, Q^j of four Q— coordinates 

in the reduced space follows from the existence of solutions for the linear (in 
Q) equation 

{Q,M} = (38) 

that is, 

dM^dQ + dM^dQ + dM^dQ = Q 
dv l dxi dx l dvi dv l dvi 
for fields without explicit time- dependence. 

For constant uniform fields one would obtain as solutions of f|39|) the 
following set of coordinates for the reduced space 



Q = x +^ [vl - u±) x B 

e Lb 

q! 0) = b-v 



(40) 
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Corrections to these coordinates for the general case arise from the devia- 
tion from uniformity in the fields. Therefore an iteration scheme may be 
devised to construct these corrections with convergence dependent on the 
fast approach to zero of the higher space derivatives of the fields. Write ( 1391) 

as 

n4~Q = o (41) 



dzi 



with 7 = \ ^ x ,j v ) and z 



Consider now an arbitrary unit vector a transverse to b- We may take 



cxb 



cxb 



c being a fixed vector in space, non collinear with b- Then define 



mc 
e\B\ 



cos 



-1 



a ■ [v x — v>x 



v± - u ± 



(42) 



/ 



6o is the angle variable conjugate to M in the case of uniform fields. Indeed, 



in this case ^i^-Oq reduces to 



v± - u± 







H B kj - 



d 



dxj mc dv, 



9 n = 1 



We now write ( )4Ti) as 







7 " (g(W + g« 



dz. 



)=0 



or 



Putting 



g (1) = -^7^q (0) 



(43) 



(44) 
(45) 
(46) 



one cancels the lowest order (in the field derivatives) terms in the right-hand 
side of ( T4"5"l) . Iterating the procedure one finally obtains 



Q 



1 + *o7i<g: 



(47) 
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As stated, convergence of the formal series (1471) will rely on the fast conver- 
gence to zero of higher order space derivatives of the fields. These reduced 
space coordinates are the variables that should enter into the gyrokinetics 
reduced Vlasov equation. 

In conclusion: By proving the existence of a vector function u±_ satisfy- 
ing Eq.( l22l we have somehow reversed the usual approach to gyrokinetics. 
Instead of the (still open) question of asymptotic gyroangle independence in 
the perturbative approach, we have shown the existence of an exact reduction 
and, whenever approximations are needed, they will focus on truncations of 
the exact u± and the corresponding approximations for the set of variables 
(Q, Qvj in the reduced space. 
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